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I - INTRODUCTION /772 » 

Thia work ia baaed on a dual hypotheaia: 

- first, reaote sensing data in thernal infrared and 
■icrowaves (surface tenperature and huaidity), such as can 
be obtained near the soil or froa an aircraft or 
satellite, appear spatially very variable at a given 
instant, on a "aeso-scale" and froa one parcel to another 
[1], but also on a aore local scale, withiu the saae 
parcel [2j; 

- second, the hydrodynaaic characteristics of soils 
(equations for pressure-huaidity h(0) and hydraulic- 
huaidity conductivity 1(8)) are well known to be very 
variable in space [3, 4]. 

Since transfers of aass and energy in soil and ataosphere 
are related, it is possible to try to discover to what extent the 
variability of hydrodynaaic characteristics is responsible for 
that of surface paraaeters. 

Next to purely experiaental approaches [5], soae tentative 
study approaches on this influence have already been used 


VNuabers in the aargin indicate pagination in the foreign text. 
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successfully with simulation models to describe the spatial 
variability of aoil characteristics using the scaling factor 
theory [7, 8]. 

The work described herein was done in this wanner, using 
previous studies on: 

- the model itself and its base equations; 

- processing results, viewed here statistically. 


II - VARIABLE HYDRODYNAMIC SURFACE CHARACTERISTICS 

The principle of conservation of mass on the soil surface 
implies algebraic equality between evaporation flux and so'l 
water flux, which are written respectively: 

e -r7T.Mp.-p.> (1) 

♦l - - * - - <K ♦ B* h ) 1C (2) 

(h v - convective exchange coefficient for water vapor, function 
of surface roughness z , wind velocity U , and the thermal 

O 8 

structure of the atmosphere; P and P - water vapor pressures on 

8 8 

the surface and at a reference atmospheric level; R y - perfect 
gas constant for water vapor; T^ - average temperature between 
these levels; z - positively counted portion of the surface; D^ 
and D _ - diffusion coefficients). 

The value of this mass flux is linked to instantaneous /773 
climatic conditions, as shown iu the energy balance equation: 

R n = L . B + H ♦ G (3) 
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R q - density of clear radiation flux; L - latent heat of water 
vaporisation; H - density of aeasurable heat flux; G - density of 
heat flux in the soil), with: 

*n ■ (l-a)Rg ♦ R, - coT g * 

H - pC p h h (T, - T.) 

c ■ • i. il . j . it 

cT iz ch 

(a and C - albedo and eaisaivity of the surface; B - solar 

g 

radiation; R - ataospheric radiation, <J - Stef an-Boltsaan 

constant; h^ - exchange coefficient for aeasurable heat; T a and 

T - surface teaperatures at the reference level; p and C 
a p 

voluaetric aass and specific heat of the air; D cT and - 

diffusion coefficients for heat). 

Thus, for cliaatic conditions (R . R , T , P , 0 ), soil 

I 8 i B fi 

characteristics (a, £, z q , K, X), and given soil conditions 
(dh/9s, 3 t/9z), surface variables (T^ and 0^ for exaaple) 
constantly assuae equilibriua values, which thus appear to be a 
function particularly of: 

- values of K and h instantaneously on the surface 

- overall behavior of the soil, which will regulate the 
aaount of water and energy on the surface and which can be 
described with the following conservation equations: 


(4) 

(5) 

( 6 ) 


<*<•> - ^ ( (K " Dyh> li * °- T li ' K) 


il . i_ { — ♦ D il 

It it it cT )s 


- /-v > T » 

il ■ il 


(7) 

( 8 ) 


( C h ( 0 > s ' capillary capacity; 0,(9) - voluaetric calorific 

capacity, D v ^* ^ c h* D cT ~ diffusion coefficients); these 

equations lead to the appearance of hydrodynaaic characteristics 
h(9) and K(6). 
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Ill - 8PATIAL VARIABILITY AND APPLICATION OF SCALING FACTOR TO 
POROUS ENVIRONMENTS 

For a given level of volumetric humidity 0, hydraulic 
conductivity K and effective pressure h vary by several orders of 
magnitude from one soil to another and even in the same soil [3]. 
If there is a certain number of study and treatment methods for 
this spatial variability [8], the approach in terns of scaling 
factor application appears to be of great interest with regard to 
the work discussed here. 

31. Introduction to the scaling factor /774 

In the case of two geometrically similar porous environments 
[9] — one environment (i) and one reference environment (A) — 
dimensional analysis of the invariability of the surface tension 
and kinematic viscosity coefficients leads to: 

h t (0) = h t (0) 

Xj" 2 K. (0) = X # ' 2 K, (0) 

( X ^ and X^ being the characteristic lengths of the internal 
geometry of the two environments). 

Thus, by assimilating environment ($) to the seeling 
average, characteristics h^(0) and K^(0) of any similar 
environment will be known if average curves h* (0) and K* (0) and 
the value of (i) in the scaling factor GT = X^/X^ are known, 
since then: 

h i (6) = *i~ l h * (0) 

K i (0) = *i K * (0) 


(9) 

( 10 ) 
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Experimentally , it is ststsd that application of the seal inf 
factor [10] , which consists of reducing variables by using the 
scaling factor, significantly reduces the dispersion of 
experimental points around the average curves for a given coil: 
a pudo logic unit could thus, in preliainary analysis, be 
considered as a group of sub-units siailar to each other (if this 
hypothesis were exact, the points would coalesce conpletely). 

The distribution of the scaling factor in the field appears 
to be noraal or, seeaingly sore often, log-normal [8]: knowledge 

of an average or a deviation-type is enough to describe it. 

32. Consequence for modeling 

The introduction of this idea of scaling factor permits 
relatively easy achievement of the objectives of this work, 
insofar as knowledge of the distribution law can supply a 
satisfactory description of the spatial variability of 
hydrodynamic characteristics. 

In this spirit, the projects already mentioned [6, 7] 
consisted of a study on sensibility of soil-atmosphere transfer 
models with variations in scaling factor, considered as an entry 
parameter: analysis was done on the pertinent exit parameters 

(temperature, humidity, evaporation). 

Here we are going to substitute a stochastic approach /775 
for this determinant approach, with a view toward supplying not a 
conservative group of results, but rather statistical 
distributions and therefore occurrence probabilities of this or 
that event; this would also make it possible to test the behavior 
of models which are determinant with regard to the stochastic 
results. 

The procedure followed will therefore be the following: 
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using a model of paired transfers of lass and energy between bare 
soil and atmosphere, sore complete than preceding models or 
■odels which are currently being experimentally validated [11], a 
satisfactory nuaber of simulations will be done for each case 
chosen (one hundred, corresponding to that aany values of the 
scaling factor taken in accordance with the distribution law) to 
be able to consider the output variables as uncertain variables 
and thu to subject then to the treatments which are standard for 
the material. 


IV - THB SIMULATION MODEL 

The aodel used simulates transfers of water (liquid or 
vapor) and energy into a soil-atmosphere system limited at its 
base by a phreatic cover and at its top by an atmospheric level 
for which the standard meteorological parameters are assumed 
known: Rg, T , P , U , and precipitation flux (duration, 

fl ® 0 

intensity). Transfers in the soil are described by two partial 
derivative equations, drawn from the theory of Philip and De 
Vries [12], formulated in pressure and temperature, and modified 
to allow for a certain number of additional phenomena [13, 14]. 
They are solved by an implicit diagram of finite differences. 

The conditions at the upper limit are supplied by simultaneous 
resolution of mass and energy balances on the surface of the 
soil. The convective exchange coefficients are formulated from 
results obtained by Dyer and Hicks [15], Webb [16], and Paulson 
[17] for stability corrections and Brutsaert [18] for roughness 
lengths. A detailed description of the model is given in [11]. 


V - SIMULATION CHARACTERISTICS 
51. Soil 
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The soil chosen for the study is PANOCHE SILT LOAN [19] 
whose average characteristics are shown in Fig. 1. Calorific 
capacity and thermal conductivity are calculated according to the 
De Vries sodel [20]. Scaling factor <X (Fig. 2) follows a log- 
normal law with parameters: 


“Log <* 


-0.616 


Log a = 116 

Albedo and enissivity, functions of soil surface /776 

humidity, are shown in Fig. 3. The length of aerodynamic 
roughness is 1 am. 


52. Climate 

Simulations are done una><sr conditions typical for the month 
of Nay in Avignon (southeaster.: France); the daily change in 
meteorological parameters is shown in Fig. 4. 

53. Initial conditions 

The initial humidity and temperature profiles are for dry 
conditions, and a 1.5-a cover is used. Since these profiles are 
not in balance with the imposed climatic conditions, two days of 
simulations are included for each case, and onl> results from the 
third day on are used. 

54 . Simulation process 

One hundred simulations, lasting one hour, are done for an 
equal number of scaling factor values in accordance with the 
distribution law. 
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VI - RESULTS 


Tbe results concern essentially surface variables 
( teaperature , huaidity) and evaporation. 

Figure 6 shows values measured at 3:00 p.m. for these 
variables as a function of the logarithm of the scaling factor. 
Low surface humidity and evaporation are immediately noted; they 
are due to the dry initial conditions. Humidity varies by more 
than 1% between the extreme values cf the scaling factor, and it 
decreases steadily as the scaling factor increases. Evaporation 
increases, but it stabilizes for high values of (t: a decrease in 

humidity following higher evaporation cauwes a decrease in 
conductivity; the scaling factor does not increase. This 
behavior is seen again for accumulated evaporation at the end of 
the day, which appears as a linear function of Log (X over a large 
part of its distribution (95*). Surface temperature decreases as 
evaporation increases; it rises again for high values of CL: 
since evaporation there is constant and huaidity decreasing, the 
energy balance can only lead to higher values. 

Therefore, evaporation end temperature distributions appear 
binodal (fig. 7). The variation coefficient is low for /777 

surface parameters (0.5* for huaidity, 4.5* for teaperature) and 
higher for evaporation (22* at 3:00 p.m., 27* for the accumulated 
value) . 

To find out whether it is possible to characterize the 
average behavior of a soil using just one value of the scaling 
factor typical of its distribution, three additional simulations 
were done with the modal (0.141), median (0.540), and average 
(1.058) scaling factor values. Figure 8 compares average 
evaporation obtained during the day with the 100 simulations to 
the evaporations obtained with these three values of d; Table 1 
shows the values, in these four cases, of accumulated 
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evaporation, aa well aa surface teaperature and huaidity at 3:00 
p.a. /. aoil whose scaling factor would be equal to the aedian of 
its diatribution appeara thua to have a behavior very near the 
average behavior of the aoil, which confiraa a previous result 
[ 21 ]. 


TABLB 1 



Accuaulated 
Kvaporation (aa) 

0 

s 

3:00 p.a. (*) 

T 

s 

3:00 p.a. (*C) 

Average 

siaulation 

0.58 

5.64 

44.66 

Deteraining 

siaulationa 




nodal 

(CL = 0.141) 

0 . 4i 

5.96 

44.92 

aedian 

(ot = 0.540) 

0.57 

5.60 

44. S7 

(<* = 1.058) 

0.67 

5.47 

44.62 


VII - CONCLUSION 

Although, for the siaulation conditions used, accuaulated 
evaporation for the day varies f'oa single to triple on the 
scaling factor variation range, teaperature and huaidity 
variations are auch saaller: 0.6*C and 1% (corresponding, 


11 



however, to e change in potential by a factor of 10). These 
deviations ere smaller than the variations coseonly observed 
under natural conditions, which seen to be linked in large part 
to surface structure) heterogeneities. It is stated elsewhere 
that the dry initial conditions give biwodal distributions of the 
exit parameters, except fox surface humidity. 

This work should Le done in other simulation /778 

conditions, particularly: 

- under humid initial conditions; 

- over a period of several days, which should allow ?or a 
much more complete approach to the relationships between 
the two variability systems. 
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Figure 1 . Average transfer characterist ica for Pancche Silt 
Loam. 



Figure 2 . Scaling factor 
distribution. 



Figure 3 . Albedo and 
eniaaivity . 



Figure 4 . Climatic 
characteristics of the 
simulations. 
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Figure 5 . Average temperature, humidity, and evaporation at 
3:00 p.m.; average accumulated evaporation for the day. 






Figure 6 . Distributions of temperature, humidity, and 
evaporation at 3:00 p.m., and of eccumulated evaporation for 
the day. 
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Figure 7 . Daily change in 
evaporation: 

- average ( ) 

- determining 

. (X modal = 0.141 ( ) 

. < t median = 0.540 
. fit average = 1.058 ( ) 
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